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A  MICROCOMPUTEP-CONTROLLEO  DIGITAL  ACOUISITIDN  AMD  ANALYSIS  SYSTEM 
FDR  THE  EXPENDABLE  BATHYTHERMOGRAPH 


bv 

Brian  Wannamaker ,  Roberto  Rossi,  Phlllo  Nesfleld  and  Paolo  Sala 


ARSTRACT 


For  a  number  of  years  the  expendable  bathythermograph  (XBT)  has  been  an 
Important  tool  for  basic  and  applied  oceanographic  survey  work.  This  has 
led  to  a  significant  body  of  work  being  done  on  the  accuracy  of  the  XBT 
system  as  a  whole.  One  of  the  most  important  limitations  has  been  the 
analogue  recording  and  manual  digitisation  of  the  profiles.  A  digital 
Interface  Is  described  that  replaces  the  standard  deck  unit  and  removes 
this  problem.  The  Interface  Is  controlled  through  an  IEEE-488  Interface 
bus  by  a  microcomputer  or  a  minicomputer.  The  data  Is  stored  with  12  bit 
precision  10  times  a  second.  As  many  as  115  profiles,  each  recording  to 
450  m  In  depth,  may  be  stored  on  a  single  cassette.  The  microcomputer  and 
plotter  may  be  used  for  a  number  of  analytical  tasks  Including  Isotherm 
following,  contour  plotting,  and  sound  speed  or  density  calculations.  The 
profiles  may  also  be  automatically  compressed  by  a  "most  reversible* 
algorithm  and  coded  to  standard  format  for  radio  transmission.  The 
bibliography  contains  presently  available  research  reports  on  the  XBT. 


INTRODUCTION 


The  expendable  bathythermograph  (XBT)  was  Initially  developed  to  quickly 
determine  the  thermal  structure  of  the  upper  ocean  from  a  ship  without 
hampering  normal  operations.  It  replaced  the  mechanical  bathythermograph 
(BT)  which,  like  present  sallnlty-temperature-depth  (STD)  and  conductivity- 
temperature-depth  (CTO)  profilers,  required  the  ship  to  heave  to  during 
lowering  or  raising  of  the  XBT.  The  resulting  thermal  profile  was  used  for 
sonar  performance  predictions. 

As  oceanographers  became  Increasingly  aware  of  the  horizontal  and  temporal 
variability  In  the  oceans  from  scattered  STD  casts,  they  realized  the  value 
of  the  XBT  for  quick,  dense  profiling  to  record  this  variability.  Thus  the 
XBT  was  used  to  fill  In  the  gaps  In  STD  surveys.  The  acoustic  community 
also  realized  the  possible  significance  of  oceanic  variability  on  sound 
propagation.  At  the  same  time  the  availability  of  Increased  computing 
power  and  sophistication  of  acoustic  models  meant  that  high  resolution 
range  dependent  data  could  be  used  If  available.  For  these  reasons  there 
was  an  Increasing  number  of  studies  of  the  accuracy  and  reliability  of  the 
XBT  system,  especially  as  the  MODE  (Mid  Ocean  Dynamics  Experiment)  and 
R0LYM00E  experiments  In  the  early  70' s. 
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These  studies  compared  BT  results  of  nearly  simultaneous  STD  and  CTD  casts. 
8ecause  the  experimenters  were  interested  in  retrieving  dynamic  heights 
from  the  data  the  Intercomparison  results  were  usually  processed  as 
differences  In  the  depth  of  the  appearance  of  specific  isotherms. 

The  BTs  themselves  have  no  pressure  sensors;  depth  is  calculated  as  a 
function  of  the  time  since  the  RT  entered  the  water.  The  coefficients  of 
this  function  are  supplied  by  the  manufacturer.  Any  error  In  temperature 
measurement  or  In  these  coefficients  will  change  the  apparent  depth  of  an 
Isotherm. 

Georgl  et  al  [1],  using  direct  calibrations  of  four  hundred  BTs,  suggested 
that  the  XBT  probe-to-probe  variability  would  be  less  than  t  0.06°C  at  the 
Q5%  confidence  level  If  the  Sippican  analogue  recorder  is  not  used.  The 
standard  deviation  Is  less  than  0.03°C  at  0°C  and  decreases  to  less  than 
0.01®C  In  the  region  of  25  to  30°C.  This  decrease  with  temperature  Is 
attributable  to  matching  of  the  individual  thermistors  at  25°C  during 
manufacture.  Thus  temperature  errors  in  the  XBTs  are  well  behaved  and  can 
be  translated  into  an  error  in  isotherm  depth  by  considering  the  local 
temperature  gradient. 

MODE  and  POLYMODE  researchers  generally  used  the  deep  T5  (1800  m)  and  T7 
(750  m)  probes.  References  [2,3,4]  supplied  evidence  of  errors  in  the 
depth  coefficients  suggested  by  the  manufacturer.  Seaver  and  Kuleshov  [5] 
were  reasonably  successful  In  modelling  this  error  and  ascribed  It  to  the 
Inverse  relationship  between  kinematic  viscosity  and  temperature  and  the 
related  changes  In  pressure  drag  on  the  falling  probe.  They  suggested  that 
weight  variations  between  probes  due  to  air  entrapped  with  the  wire  should 
result  In  depth  errors  with  a  standard  deviation  of  only  1 .6  m  at  750m.  On 
the  other  hand  changes  in  surface  roughness  on  the  BT  nose  section 
sufficient  to  vary  the  point  of  transition  of  laminar  turbulent  flow  from 
the  nose  of  the  probe  to  the  point  of  maximum  cross-sectional  area  would 
create  a  standard  deviation  in  apparent  isotherm  depth  of  6.2  m  at  750  m. 

A  statistical  analysis  of  nearly  two  thousand  T4  XBTs  launched  and  recorded 
by*  O.S.  Navy  ships  has  been  reported  by  Anderson  [6].  Only  80*  of  the 
temperature  profiles  from  these  casts  were  visually  acceptable  and  40*  of  a 
Subset' Of  518  casts  met  the  accuracy  criteria  he  set.  In  commenting  on 
Anderson's  report,  the  manufacturer,  Sippican,  suggested  that  many  of  the 
failures  may  have  been  due  to  aging  of  the  probes  beyond  the  two  year 
warranty  period.  This  period  Is  set  by  the  time  needed  to  degrade  the 
Insulation  on  the  wire  and  thermistor.  A  Netherlands  Navy  study  reported 
that  the  signal  wire  becomes  sticky  because  of  age  and  high  temperatures. 
This  stiakloess  leads  to  Increased  tension  wire  stretch,  depth  errors  and 
spikes  on  the  temperature  depth  record  during  deployment  (report  of  the  9th 
MJtOC  Meeting) .  Sippican  states  that  they  will  provide  replacement  units 
If  there  Is  an  abnormally  high  failure  rate  beyond  10%.  A  review  of  the 
errors  commonly  found  with  the  standard  XBT  system  is  given  by  Kronen  and 
ftluimthal  [7], 

It  became  clear  to  many  researchers  that  much  of  the  error  In  XBT  records 
•ere-  the  result  of  the  stendard  deck  unit  and  analogue  recorder  and  the 
aemual  digitisation.  A  study  of  two  US  Navy  exercises  showed  that  more 
than  60*  of  the  BTs  collected  were  Improperly  encoded  [8],  This  deficiency 
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could  be  avoided  by  digitally  recording  the  data.  The  manufacturer  has  for 
some  years  offered  a  digital  shaft  encoder  or  retransmitting  potentiometer 
mounted  on  the  analogue  recorder  temperature  servoassembly.  A  system  at 
SACLANTCEN  used  the  encoder  for  its  first  digital  interface  [9].  This  was 
successful  In  alleviating  time  consuming  and  error  prone  pre-analysis 
procedures  but  was  limited  In  precision  to  0.05®  because  of  the  encoding 
system.  The  potentiometer  was  used  as  part  of  a  system  by  McDowell  and 
Oorson  [3]  along  with  a  12  bit  analogue-to-dlgital  converter  and  cassette 
recorder.  This  was  the  basis  of  a  commercial  system  by  Bathy  Systems  Inc. 

Both  systems  were  constrained  by  the  limits  of  the  signal  processing  before 
the  digitisation.  Other  systems  were  built  to  completely  replace  the 
standard  on-deck  electronics  (e.g.  [10])  and  some  were  controlled  by 
microprocessors  [11].  In  1985  Slpplcan  has  also  announced  a  digital 
Interface  (MK-9)  with  attached  microcomputer  (HP-85)  although  this  system 
records  only  one  XBT  profile  per  cassette. 

A  specific  system  to  collect  temperature  gradient  measurements  along  with 
temperature  measurements  has  been  developed  by  Lange  and  Johnson  [12] .  It 
replaces  the  XBT  with  a  recoverable  instrument  package  for  rapid  profiling 
that  Is  being  developed  at  the  University  of  Washington  [13]. 

This  memorandum  discusses  a  microcomputer  controlled  deck  unit  system 
developed  at  SACLANTCEN.  As  yet  It  has  not  been  compared  to  other  systems 
of  the  same  concept.  From  the  available  literature  It  appears  that  the 
SACLANTCEN  system  offers  greater  data  storage  than  some,  easier  one-person 
use  and  more  complete  Internal  self-checking  and  flexibility  than  others. 

The  first  section  of  this  report  will  outline  the  use  of  this  system  as  a 
complete  unit.  This  will  be  followed  by  more  detail  about  the  hardware  and 
software.  References  concerning  XBTs  have  been  extracted  from  a 
computerized  bibliography  accessible  on  the  SACLANTCEN  main  frame  computer 
and  Is  provided  at  the  end  of  this  memorandum. 


1  OPERATION  AS  A  SYSTEM 

The  components  of  the  XBT  acquisition  and  analysis  system  are  shown  In 
Fig.  1. 

The  controller  Interface  (C/I)  completely  replaces  the  normal  deck  unit; 
the  cable  from  the  launcher  to  the  normal  unit  Is  simply  disconnected  and 
reconnected  to  the  new  controller  Interface  (C/I)  by  the  same  type  of 
amphenol  connector.  (This  Is  especially  Important  to  take  advantage  of 
ships  of  opportunity.  One  of  the  units  has  already  been  used  aboard  the 
Italian  Navy  hydrographic  ship  MAGNAGHI  (IT)  In  1984.)  Hand  launchers  may 
be  connected  via  a  rear  terminal  strip.  The  C/I  Is  In  turn  controlled  by  a 
micro-  or  mini-computer. 

One  person  can  use  the  system  If  necessary;  communication  with  an  assistant 
at  the  launcher  Is  not  needed.  When  the  acquisition  program  runs,  the 
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OPTIONAL  FLOPPY  OISC 


FIG.  1  EXPENDABLE  BATHYTHERMOGRAPH  (XBT) 
ACQUISITION  AND  ANALYSIS  SYSTEM 


operator,  when  asked,  reads  the  real  time  to  form  part  of  the  header  record 
of  each  data  file. 

The  labels  of  the  eight  special  function  keys  are  displayed  along  the 
bottom  of  the  computer  monitor  as  shown  in  Fig.  2.  These  keys,  always 
active,  are  programmable  to  transfer  program  execution  to  numbered 
subroutines.  For  example  "UNIT"  initialises  the  output  plot  parameters. 
The  operator  is  also  asked  for  the  BT  type  since  the  fall  rate  of  the  probe 
and  data  file  size  required  are  type-dependent.  These  values  are  then 
available  for  use  In  the  routines  until  the  program  is  halted.  When  the 
parameters  have  been  Initialized  the  operator  presses  “GO"  and  the  CRT 
prompts  "PLEASE  RELOAD"  (the  XBT  launcher). 

If  working  alone  the  operator  can  now  go  to  the  launcher  and  load  and 
launch  the  BT;  Internal  checks  and  data  measurement  are  automatic.  The 
system  could  be  set  In  a  loop  to  prepare  immediately  for  another  BT  while 
the  operator  remained  at  the  launcher.  The  data  is  plotted  on  the  computer 
monitor  during  acquisition.  The  conversion  from  probe-thermistor 
resistance  to  temperature  is  through  the  polynomial. 


***  **  **.* 

*  *  *  t  * 

*  *  *  * 

*  *  85  *  *** 

X  t  *.  * 

*  *  *  *  .* 

***  it  *** 


TEST  AG  PK  WESG 

IN  I T  GO  FLOT  REPLAY 


FIG.  2  COMPUTER  MONITOR  SHOWING  THE 
EIGHT  SPECIAL  FUNCTION  KEYS 
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T  =  63.2204  -  11.9493R  +  1.10525  R2  -  0.0552171  R3  +  0.0010914  R4 

which  was  fitted  to  the  calibration  data  supplied  by  the  manufacturer. 

Temperature  is  in  degrees  Celsius  (BC)  and  resistance  is  in  kfl.  The  depth 
calculation  is  done  with  the  formula: 

D  =  (6.472  -  0.00216  *  J/10)*  J/10 

for  T4,  T7  and  T2  type  XBT's  and 

D  =  (6.828  -  0.00182  *  J/10)  *  J/10 

for  T5  uni  ts 

with  D  the  depth  in  metres  and  J  the  sample  number  (1  sample  is  equivalent 
to  0.1  s)  since  entering  the  sea. 

Wire  breakage  during  the  acquisition  is  reported  by  the  C/I  status  word. 

The  data  from  forty-two  T4  BTs  will  fit  on  one  DC  100A  data  cassette.  This 
is  a  computer  operating  system  limit  and  could  increase  with  software 
changes  to  115  profiles. 


2  HARDWARE 

The  Controller  Interface  (C/I)  replaces  the  standard  Sippican  desk  unit  and 
is  transportable,  weighing  only  about  9  kg.  The  case  can  be  mounted  In  a 
rack. 

The  standard  deck  launcher  or  hand  launcher  may  be  connected  either  via  a 
8  pin  amphenol  connector  or  a  terminal  strip. 

The  microcomputer  and  XBT  C/I  are  linked  by  an  IEEE-488  or  GP-IB  (General 
Purpose  Interface  Bus)  also  referred  to  as  HP-IB.  A  wide  range  of  micro 
and  mini  computers  can  thus  be  used.  However,  SACLANTCEN  uses  an  HP -85 
micro-computer  or  HP  21MX  series  mini -computer. 

The  computer  sends  control  words  to  the  C/I  via  this  link  and  can  request 
status  reports  from  It  as  well  as  data.  The  C/I  also  displays  its  own 
status  and  that  of  the  HP-IB  as  well  as  the  probe  type  in  use  via  front 
panel  light  emitting  diodes  (LEDs). 

A  simplified  block  diagram  of  the  C/I  circuit  is  shown  in  Fig.  3. 


2.1  The  Self  Test  Circuit 

This  circuit  checks  the  launcher,  probe  and  laun-:her-to-C/I  cable  each  time 
the  breach  of  the  launcher  Is  opened  and  closed. 


The  circuit  tests  for  certain  voltage  thresholds  above  which  the  probe  and 
connections  are  indicated  to  be  good.  If  this  test  is  completed  the  input 
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FIG.  3  BLOCK  DIAGRAM  OF  THE  XBT  CONTROLLER  INTERFACE 


terminals  of  the  C/I  are  connected  across  a  10  kfl  precision  resistor.  Data 
are  sent  to  the  host  for  1  second,  enough  time  for  100  samples, and  are 
checked  for  a  temperature  of  9.92  ±  0.01“C.  (This  temperature  Indicates 
that  the  analogue  interface  is  working  properly.  Any  difference  from  this 
'tlue  represents  a  bias  that  is  correctable  later  during  the  data 
analysis).  The  host  computer  can  then  send  a  control  word  to  enable  the 
launch. 


2.2  Resistance  to  Voltage  Converter 

After  the  completion  of  the  check  and  the  receipt  of  an  "ENABLE"  command 
the  probe  is  connected  to  the  input  of  the  resistance-to-voltage  converter 
shown  in  Fig.  4. 

This  portion  of  the  C/I  is  as  described  in  Stegen  et  al  [10].  The  probe 
thermistor  forms  one  leg  of  a  self-balancing  bridge  circuit.  The  bridge  Is 
formed  by  the  thermistor,  two  precision  resistors  and  J-FET  Tl.  The 
resistance  of  Tl  is  constrained  to  equal  the  thermistor  resistance  by  the 
feed  back  of  operational  amplifier  Al.  Differential  amplifier  A2  provides 
an  output  voltage  directly  proportional  to  the  resistance  of  Tl  and  thus  of 
the  probe  thermistor. 
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. .1  Analogue  to  Digital  Conversion 

The  analogue*  ou of  A?  is  fed  to  the  Input,  of  a  1A  hit  analogue  to 
coital  converter,  the  most  significant.  ]?  hits  abut  the  output  of  A 
nyclir:  a  nit  conn  if”,  tboy  cron  ting  the  10  bit  data  word  that  in  sent  to  the 

I'omiiutnr  1  hr  .i.i r.»  :  .rivet  r.1  on  in  completed  in  100  ms  to  an  accuracy  of 

"ht  coimuM'  : s  v<i  to  toro  h>/  the  closing  of  the  sensor  circuit  by 

sea  |'i 1  cn  when  ':••((’  probe  enters  the  water.  This  cou ate;*  is  use!  as  a 
dingn-'-.tSs  tool  U;  df-frrmi ne  that,  nn  data  samples  are  helm  lost  between 


■  if  f;  ;><  ‘  and  Control  Logic 

'“ns  He-iii  piuKessoi  a  Single  LSI  monolithic  chip  handl Inn  all  the  data 
o 1 oo  in  both  di ree  tiers  h? tween  the  host  computer  and  the  XOT  C/I  using  the 
standard  n-.i.h-dUfl  i i > )  protocol .  The  associated  control  logic  Insures 

that  all  i'; .  id'  e.,!.  1  m  signals  'Yom  and  to  the  lib- IP,  processor  are  handled 

prop/,  pi  y  . 


Or,  r,ardpe  !y  e nonce ; 

fho  samp  in  vguonce.  supplies  the  analogue  to  digital  trigger  signals  at  a 
frequency  •  i  y  Si/,  fieri ved  by  electronically  dividing  the  output  Of  a 
precision  1  •"'/  •".•yv.tni  oscillator.  This  circuit  also  controls  the 
genera t< or  ni  tug  el  t  test  digital  data  and  the  four  hit  counter, 
Eventually  f  the  i.  will  incorporate  the  ah  ill  tv  to  change  sampling  rate 

coder  ho  ■;  |  r  11 ' 


Test  pci  tor* 

The  test  tea  Mir*’  allows  a  check  to  be  node  of  the  C/1  itself  and  also 
gonj^ratos  fi':  >  chirking  of  the  computed  acquisition  and  processing. 

T  iv-re  are  fb-i  g  1 1  •  oi'ent  'rgfs: 

•'  Profile  simulation:  this  test  generates  256  samples  that  simulate  a 

r us  1 1. 1: ah f .1 :  c iiTiT i g e  (ranging  from  VM.l  to  lR090.fl  ohms  or  !M.?11*C  to 
1  .  The  resulting  profile  is  shown  in  fig.  5. 

M  Minimum  re-sl stance  1  n«- :  this  test  nrodiif.es  an  output  of  T?f>1.7  ohms 


Maximum  rcsisi  an  re  value :  this  test  produces  an  output  of  lfl090.11  n 

•j’  y  a l ’•*  /'"  *1 

In  each  ease  Um  te  .1  i‘  programmed  by  the  host  computer  and  Initiated  by 

'k'i't  ’  eg  fhi‘  j;r‘  k  boft.oyi  on  the  front  panel. 


BEST  AVAILABLE  COPY 
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FIG.  5 

SIMULATED  TEST 
PROFILE 


2.7  Accuracy 

Figure  6  is  a  plot  of  the  output  resistance  value  of  the  controller/ 
interface  measured  by  the  micro-computer  vs  the  voltage  output  of  amplifier 
A2  (see  Fig.  4).  This  line  is  linear  with  a  variation  equal  to  ±  1  bit  of 
the  analogue/digital  (A/D)  converter. 

Figure  7  is  the  result  of  an  attempt  to  measure  the  absolute  accuracy  of 
the  C/I.  The  XBT  probe  was  replaced  with  a  resistor  decade  box  (which  has 
a  precision  of  0.5%).  The  difference  between  the  output  resistance  value 
and  the  decade  value  is  plotted  over  the  full  range  of  the  output  of 
amplifier  A2.  Howaver,  at  the  bottom  of  the  scale  the  ’error''  is  less  than 
the  uncertainty  in  the  input  value. 


FIG.  6 

OUTPUT  RESISTANCE 
OF  CONTROLLER 
INTERFACE  VS 
VOLTAGE  OUTPUT  OF 
MEASUREMENT 
AMPLIFIER 
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FIG.  7  XBT  RESISTANCE  MEASURING  ERROR 

The  absolute  accuracy  of  the  controller  inter¬ 
face  is  greater  than  that  of  the  standard 


2.8  Navigation  Satellite  Interface 

The  SACLANTCEN  research  ship  MARIA  PAOLINA  G.  is  equipped  for  satellite  and 
Loran-C  navigation.  The  computer  controlling  this  can  be  linked  to  the 
HP-85  through  an  Interface  built  at  the  Centre.  Thus  the  accurate  time  and 
geographical  coordinates  can  be  included  in  each  XBT  data  file.  This 
Interface  will  be  discussed  elsewhere  in  connection  with  a  microcomputer 
system  that  collects  meteorological  data. 


3  SOFTWARE 

The  flexibility  of  the  HP-IB  interface  bus  and  of  the  HP-85  microcomputer 
has  made  the  programming  of  the  data  acquisition  relatively  easy. 
Acquisition  routines  for  the  HP21-MX  mini  computer  series  have  also  been 
written  to  Incorporate  the  XBT  data  into  the  Integrated  Oceanographic 
System  (I0S).  The  microcomputer  allows  a  more  general  use  of  the 
advantages  of  the  new  interface  since  it  is  Independent  of  the  110V 
generators  and  other  computer  processing  aboard  ship. 

The  basic  HP85  includes  16Kbytes  (14  576  are  available)  of  read/write 
memory  with  an  Internal  operating  system  running  on  an  expanded  version  of 
ANSI  BASIC,  an  integral  cassette  drive,  a  CRT  screen  and  a  thermal  printer. 
The  HP85  used  In  the  XBT  system  has  been  expanded  with  plug-in  units  for  an 
additional  16K  of  memory  and  for  external  printing  and  plotting. 
Input/output,  matrix  functions  and  mass  storage.  However  the  matrix  and 
mass  storage  functions  are  not  presently  required  for  the  data  acquisition 
or  analysis. 

Two  separate  routines  have  been  developed;  the  first  Is  used  for  data 
acquisition  and  single  profiles  and  the  second  is  concerned  with  stored 
data. 

The  first,  besides  acquiring  the  profile,  can  run  Interface  checks,  plot 
the  profile  (as  acquired  or  from  tape),  compress  the  data  and  print  the 
profile  in  a  standard  format,  as  determined  by  the  operator. 
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Flags  (variables  with  value  1  or  0)  control  optional  program  branching  that 
is  generally  held  constant  for  a  run.  This  includes  such  things  as 
automatic  hardcopy  of  profile,  type  of  surface  corrections  etc.  These 
values  are  read  from  data  statements;  in  BASIC  they  can  be  quickly  changed 
by  the  operator. 

Because  of  the  small  size  of  the  computer,  efforts  were  made  to  save 
storage  space  which  was  shared  by  proqram  and  data.  The  memory 
requirements  for  array  variables  in  the  HPR5  are  shown  below  in  Table  1. 


TARLE  1 

Memory  Requirements,  HP-85 


Type  of  Array  Variable 

Bytes  of  Memory 

Full  Precision 

Short  Precision 

Integer 

String 

8  +  8/element 

8  +  4/element 

8  +  3/element 

8  +  1/element 

Depth  information  is  not  stored.  It  depends  only  on  BT  type  and  sample 
number.  The  BT  type  is  kept  as  part  of  the  file  name  and  used  to  branch  to 
the  appropriate  depth  calculation  during  analysis.  Although  not  a  problem 
so  far,  loss  of  data  by  deficiencies  in  the  system  can  be  determined  by 
following  the  consecutive  counter.  The  data  arrive  at  the  microcomputer  as 
16  bit  Integer  values  that  represent  an  unsealed  resistance.  Without  any 
loss  of  accuracy  this  can  be  converted  to  a  pair  of  characters,  8  bits  at  a 
time. 

Tape  files  are  created  just  before  data  acquisition  because  there  is  no 
quick  way  to  determine  if  an  existing  data  file  is  full  or  empty.  Storing 
resistance  values  rather  than  calculated  temperatures  is  consistent  with 
the  policy  of  the  SACIANTCEN  Applied  Oceanographic  Group  to  keep  raw  data, 
or  data  that  can  be  unambiguously  reconverted  to  that  form.  In  this  case 
the  resistance  to  temperature  conversion  is  non-linear.  However,  if  this 
relationship  was  later  found  to  be  in  error  or  if  an  improved  expression 
was  determined,  temperatures  already  converted  could  not  be  corrected 
whereas  the  raw  data  would  still  be  available  for  conversion. 

Once  the  program  is  running  the  HP-85  becomes  a  dedicated  computer.  To 
prevent  Inadvertent  program  halts  the  keyboard  becomes  inactive  except  for 
certain  specific  keys. 

The  flexibility  of  the  software  control  and  HP-IB  allows  the  correction  of 
other  possible  errors.  These  errors  include  initiating  acquisition  with 
"GO"  twice  or  opening  the  launcher  breach  between  loading  and  pulling  the 
releasing  pin  of  the  probe. 
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The  data  transfer  Is  done  In  pieces  and  interrupts  other  executing  tasks. 
However  no  data  are  lost  if  incoming  information  is  converted  and  plotted 
as  a  temperature  vs  depth  profile  during  acquisition.  The  Internal  screen 
has  a  horizontal /vertical  resolution  of  256  «  182  dots  and  the  data  are 
scaled  between  the  limits  set  during  the  initialisation  (INIT)  phase.  A 
portion  of  a  profile  may  be  shown  at  full  resolution  or  a  piece  of  it  may 
be  displayed  on  an  external  plotter,  also  connected  throuqh  an  HP-IB  link. 

A  T-4  BT  operating  for  90  s  supplies  about  900  temperature  data  points. 
Many  of  these  data  are  redundant,  especially  those  in  the  surface  mixed 
layer  and  in  the  generally  isothermal  deep  layers.  The  SACLANTCEN  data 
base  has  a  system  limit  of  125  temperature  depth  pairs;  ships  reporting 
temperature  profiles  to  forecast  centres  by  radio  are  limited  to  20  data 
pairs  in  a  specified  format.  Obviously  some  way  to  compress  the  data  is 
necessary. 

This  may  be  done  by  someone  reading  selected  values  from  the  trace  or 
through  the  use  of  a  digitising  table,  both  procedures  are  tedious  and  thus 
error-prone.  Also  the  digiting  must  properly  take  into  account  the  non¬ 
linearity  of  the  chart  scaling.  In  any  case  the  omanufacturer  attributes 
half  of  the  limit  of  the  system  uncertainty  of  0.2  to  the  analogue  chart. 

The  system  described  here  uses  a  method  developed  by  R.  Wlnterburn  and 
Wannamaker  [14].  This  algorithm  aims  to  determine  a  reduced  set  of 
temperature,  depth  pairs  that  when  joined  with  straight  lines  will 
reconstruct  the  input  profile  within  a  known  limit. 

Basically  the  algorithm  chooses  as  significant  points  those  at  which  the 
change  in  temperature  gradient  is  above  a  set  threshold.  If  so  the  point 
is  kept  as  a  significant  point  for  the  reconstruction  of  the  profile;  if 
not  it  is  subjected  to  the  following  test:  the  true  value  at  the  midpoint 
between  this  point  and  the  last  significant  point  is  checked  to  determine 
if  it  is  within  a  preset  limit  of  the  straight  line  joining  these 
significant  points.  If  not,  the  midpoint  becomes  significant  as  well.  The 
international  standard  for  the  exchange  of  oceanoqraphic  data  IOC  [15] 
states  that  for  XBT  data: 

...  flexure  points  (should  be)  determined  in  such  a  way  that  linear 
Interpretations  fall  within  ±  0.2°C  of  the  original  record" 

Wlnterburn  and  Wannamaker  [14]  discussed  this  method  for  200  profiles  from 
the  Alboran  Sea  and  Gulf  of  Cadiz  in  frontal  regions.  In  the  Gulf  of 
Cadiz  region  the  profiles  show  a  sudden  increase  in  temperature  with  depth 
when  the  probe  entered  the  Mediterranean  outflow.  He  found  that  a 
threshold  value  of  0.035°C/m  for  the  gradient  change  reduced  the  original 
number  of  points  of  the  profiles  by  85%.  A  check  against  the  original  data 
at  each  original  depth  value  showed  an  average  standard  deviation  of 
0.022#C  and  maximum  error  of  0 . 17°C . 

Although  a  threshold  value  of  0.3°  c/m  was  required  to  reduce  each  of  15 
Alboran  Sea  profiles  to  less  than  20  points  on  the  HP  85,  the  error,  when 
compared  with  the  original  profile,  was  generally  less  than  0.2°C. 

Surface  temperature  measurements  present  another  aspect  of  the  XBT. 
Results  calculated  from  modelling  the  time  response  of  the  sensor  have 
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shown  that  an  XBT  cannot  give  an  accurate  temperature  reading  until  about 
4  m  depth  [16].  Stegen  et  al  [10]  suggest  Ignorlnq  the  data  during  the 
first  four  time  constants.  However,  some  researchers  have  found  good 
comparisons  with  simultaneous  measurements  by  other  Instruments  T 1 7 1 .  In 
the  present  system  the  user  may  opt  through  Hag  setting  to  Ignore  the 
upper  3.7  m  reported  by  the  XRt  or  to  extrapolate  back  up  to  the  surface 
using  the  slope  between  6  and  4  m  [16’. 

The  resulting  reduced  profile  may  then  be  automatically  coded  into  standard 
radio  message  format  described  by  Fisher  [18].  The  depth  values  are 
rounded  off  to  the  nearest  metre  and  the  temperature  to  tenths  of  a  degree 
as  part  of  the  coding. 

An  example  profile,  packed  profile  and  message  are  shown  in  Fig.  8.  At 
least  one  other  compression  algorithm  is  in  use  for  XBT  data  [see  Mesecar 
and  Wagner,  [11].  Other  data  compression  schemes  exist  and  research  is 
continuing  to  determine  the  one  most  suitable  for  this  application. 

A  reduced  data  base  of  T-S  diagrams,  although  not  implemented  in  the 
system,  Is  suggested  for  future  development,  ilean,  median  or  some  other 
representation  of  the  T-S  relationship  in  a  geographical  area  at  a  given 
time  of  year  could  be  stored  on  a  floppy  disk  with  the  appropriate  hardware 
expansion.  Flier!  and  Robinson  [19]  have  suggested  an  interpolation  scheme 
with  a  number  of  T-S  curves  rather  than  one  representatl ve  curve.  Over 
limited  regions  or  within  limited  periods  or  if  polynomial  approximations 
for  the  relations  were  used  the  T-S  curves  might  fit  In  the  computer  memory 
along  with  the  program.  These  data  bases  would  allow  more  accurate 
estimates  of  the  salinity  profile  and  thus  density,  sound  velocity  and 
geostrophic  currents.  A  data  set  of  mean  T-S  diagrams  has  been  developed 
for  5°  latitude/longitude  squares  for  the  M.  Atlantic  and  Pacific  [70]. 
One  for  the  Medl terranean  could  be  determined  from  the  SMODS  data  base  at 
SACLANTCEN  [21].  A  non-regular  grid  is  suggested  so  that  separate  files 
would  exist  for  Inside/outside  the  Alboran  Sea  gyre  for  example.  Such 
processing  has  been  done  with  BT  data  but  usually  only  after  the  cruise. 
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Once  a  series  of  XBT  profiles  are  on  tape  they  may  be  processed  quickly  in 
sequence.  Assuming  they  are  all  from  the  same  BT  type  a  depth  array  can  be 
filled  once  and  retained  for  use  in  the  second  basic  software  routine. 
Corresponding  temperatures  are  calculated  for  each  data  file  upon  reading 
and  stored  as  short  precision  values  (See  Table  1). 

One  type  of  sequential  processing  is  sequential  plotting  (see  Fig.  9). 
Another  is  contour  plotting;  however,  this  process  may  be  too  complex  to  be 
efficient  on  a  machine  of  this  size.  The  basic  problem  is  the  decision 
process  required  to  join  contours  through  temperature  inversions  where  the 
same  Isotherm  occurs  more  than  once  on  some  but  not  all  sequential 
profiles.  For  the  present  system  this  is  left  to  human  decision.  The 
system  plots  the  depth  of  each  isotherm  with  specific  symbols.  The  user 
need  enter  only  the  contour  interval,  D.  Isotherm  determination  is  made 
easier  with  the  Intrinsic  integer,  sign  and  logical  evaluation  routines 
Inherent  in  the  HP85.  If  two  sequential  data  points  are  (Tl,  PI)  and 
(T2,  P2)  then  the  number  of  isotherms  between  them  is 

INT  (II)  -  INT  (li)  (Eq.  4) 

0  0 


13  l^>  If  19  ?1  Temperature  (’C) 


100  •  1  1  1  4  300 


ri('.  *  <Er>ri:\’TiAi.  i  lots  >t  xht  won  iff 
(')ffs-r  t".  .'°c') 


The  Isotherms  can  be  determined  by  linear  interpolation  through  a  loop. 

The  first  Isotherm  of  the  loop  is 

INT  ((Tl  -  Bl/Cl  ♦  f  f  T2  -  Tl )>  0)  *C1) 

The  underlined  expression  Is  1  when  true  and  0  when  false.  B1  Is  a  bias 
used  to  reduce  the  number  of  symbols  required.  For  the  Mediterranean  it  Is 
set  equal  to  13°C.  The  direction  of  the  loop  is  set  by  the  sign  of 
(T2-T1). 

Each  Isotherm  has  Its  own  distinctive  symbol  that  is  stored  in  a  character 
string  and  can  be  treated  as  a  one-dimensional  array.  The  correct  symbol 
can  then  be  retrieved  for  plotting  by  determining  the  appropriate  Index. 
The  symbols  used  are  limited  to  those  in  the  character  set  Installed  In  the 
external  plotter.  A  sample  plot  is  shown  In  Fig.  10a,  with  the  manually 
contoured  version  In  Fig.  10b.  Obviously  the  software  described  Is  only  a 
small  subset  of  what  could  be  Implemented  depending  upon  the  research  need. 
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CONTOUR  INTERVAL  «  1°  C 


200  J 

FIG.  10 


b) 


SAMPLE  PLOTS  OF  TEMPERATURE  PROFILES 

a)  Original  plot  with  symbolic  delineation  of  isotherms 

b)  Manually  contoured  isotherms {from  Fig.  10a). 

The  scale  at  the  top  represents  the  numbering  of  the  XBT's 
that  were  deployed  at  regular  intervals . 
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SUMMARY 


The  SACLANTCEN-developed  flexible  XRT  diqital  data  acquisition  and  analysis 
system  has  been  described.  It  is  transportable,  thus  suitable  for  ships  of 
opportunity  as  well  as  research  vessels.  As  presently  configured  it  is 
microcomputer-control  led  but  can  be  connected  to  any  computer  that  will 
support  the  IEEE-488  General  Purpose  Interface  Rus.  It  acquires  the  data 
from  the  probe  with  12  bit  precision  at  in  Hz.  The  data  can  be 
automatically  reduced  to  20  or  fewer  temperature  depth  pairs  and  coded  in 
standard  format  for  radio  messages.  The  same  concept  could  be  used  as  an 
Interface  for  the  AX8T  and  newer  XSV.  A  discussion  of  other  work  on  the 
XBT  as  an  instrument  is  included  and  a  biblioqraphy  of  pertinent  reports  is 
gi ven . 
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